in 



Search for the b — ► dj process 

D. Mohapatra, 47 M. Nakao, 7 B. D. Yabsley, 47 K. Abe, 7 K. Abe, 40 H. Aihara, 42 K. Arinstein, 1 Y. Asano, 46 
T. Aushev, 11 T. Aziz, 38 S. Bahinipati, 4 A. M. Bakich, 37 E. Barberio, 18 M. Barbero, 6 I. Bedny, 1 U. Bitenc, 12 
I. Bizjak, 12 A. Bondar, 1 A. Bozek, 24 M. Bracko, 7 - 49 - 12 J. Brodzicka, 24 T. E. Browder, 6 Y. Chao, 23 A. Chen, 21 
W. T. Chen, 21 B. G. Cheon, 3 R. Chistov, 11 S.-K. Choi, 5 Y. Choi, 36 A. Chuvikov, 32 J. Dalseno, 18 M. Danilov, 11 
M. Dash, 47 J. Dragic, 7 S. Eidelman, 1 T. Gershon, 7 G. Gokhroo, 38 B. Golob, 17 ' 12 A. Gorisek, 12 T. Hara, 29 
K. Hayasaka, 19 M. Hazumi, 7 L. Hinz, 16 T. Hokuue, 19 Y. Hoshi, 40 S. Hou, 21 W.-S. Hou, 23 Y. B. Hsiung, 23 
T. Iijima, 19 K. Ikado, 19 A. Imoto, 20 K. Inami, 19 A. Ishikawa, 7 M. Iwasaki, 42 Y. Iwasaki, 7 J. H. Kang, 48 J. S. Kang, 14 
H. Kawai, 2 T. Kawasaki, 26 H. R. Khan, 43 H. J. Kim, 15 S. M. Kim, 36 K. Kinoshita, 4 P. Krizan, 17 - 12 P. Krokovny, 1 
R. Kulasiri, 4 C. C. Kuo, 21 A. Kuzmin, 1 Y.-J. Kwon, 48 G. Leder, 10 S. E. Lee, 34 T. Lesiak, 24 J. Li, 33 A. Limosani, 7 
S.-W. Lin, 23 D. Liventsev, 11 G. Majumder, 38 F. Mandl, 10 A. Matyja, 24 Y. Mikami, 41 W. Mitaroff, 10 H. Miyata, 26 

G. R. Moloney, 18 T. Nagamine, 41 Y. Nagasaka, 8 E. Nakano, 28 H. Nakazawa, 7 S. Nishida, 7 O. Nitoh, 45 T. Nozaki, 7 
S. Ogawa, 39 T. Ohshima, 19 S. Okuno, 13 S. L. Olsen, 6 Y. Onuki, 26 W. Ostrowicz, 24 H. Ozaki, 7 P. Pakhlov, 11 

H. Park, 15 R. Pestotnik, 12 L. E. Piilonen, 47 Y. Sakai, 7 T. R. Sarangi, 7 N. Sato, 19 N. Satoyama, 35 T. Schietinger, 16 
O. Schneider, 16 J. Schumann, 23 A. J. Schwartz, 4 K. Senyo, 19 M. E. Sevior, 18 H. Shibuya, 39 A. Somov, 4 N. Soni, 30 

R. Stamen, 7 S. Stanic, 27 M. Staric, 12 T. Sumiyoshi, 44 F. Takasaki, 7 K. Tamai, 7 N. Tamura, 26 M. Tanaka, 7 
G. N. Taylor, 18 Y. Teramoto, 28 X. C. Tian, 31 K. Trabelsi, 6 S. Uehara, 7 T. Uglov, 11 K. Ueno, 23 S. Uno, 7 
G. Varner, 6 S. Villa, 16 C. C. Wang, 23 C. H. Wang, 22 Y. Watanabe, 43 Q. L. Xie, 9 Y. Yamashita, 25 M. Yamauchi, 7 
Heyoung Yang, 34 J. Ying, 31 C. C. Zhang, 9 J. Zhang, 7 L. M. Zhang, 33 Z. P. Zhang, 33 V. Zhilich, 1 and D. Ziircher 16 
i (The Belle Collaboration) 

^ , 1 Budker Institute of Nuclear Physics, Novosibirsk 

C ■ 2 Chiba University, Chiba 

Chonnam National University, Kwangju 
y ' University of Cincinnati, Cincinnati, Ohio 45221 
5 Gyeongsang National University, Chinju 
6 University of Hawaii, Honolulu, Hawaii 96822 
7 High Energy Accelerator Research Organization (KEK), Tsukuba 
8 Hiroshima Institute of Technology, Hiroshima 
^ ' 9 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 

Institute of High Energy Physics, Vienna 
• 11 Institute for Theoretical and Experimental Physics, Moscow 

| 12 J. Stefan Institute, Ljubljana 

, 13 Kanagawa University, Yokohama 

14 Korea University, Seoul 
15 Kyungpook National University, Taegu 
. 16 Swiss Federal Institute of Technology of Lausanne, EPFL, Lausanne 

5^ ' 17 University of Ljubljana, Ljubljana 

University of Melbourne, Victoria 
19 Nagoya University, Nagoya 
Nara Women's University, Nara 
21 National Central University, Chung-li 
22 National United University, Miao Li 
23 Department of Physics, National Taiwan University, Taipei 
2J *H. Niewodniczanski Institute of Nuclear Physics, Krakow 
25 Nippon Dental University, Niigata 
26 Niigata University, Niigata 
27 Nova Gorica Polytechnic, Nova Gorica 
28 Osaka City University, Osaka 

29 Osaka University, Osaka 
3C ' Pan jab University, Chandigarh 
31 Peking University, Beijing 
32 Princeton University, Princeton, New Jersey 08544 
33 University of Science and Technology of China, Hefei 
3Jt Seoul National University, Seoul 

35 Shinshu University, Nagano 
36 Sungkyunkwan University, Suwon 
37 University of Sydney, Sydney NSW 



o 
in 
o 



2 



Tata Institute of Fundamental Research, Bombay 
36 ' Toho University, Funabashi 
^"Tohoku Gakuin University, Tagajo 
* Tohoku University, Sendai 
■* Department of Physics, University of Tokyo, Tokyo 
* 'Tokyo Institute of Technology, Tokyo 
^ Tokyo Metropolitan University, Tokyo 
Jft 'Tokyo University of Agriculture and Technology, Tokyo 
J * 6 University of Tsukuba, Tsukuba 
"* Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24-061 

^ 8 Yonsei University, Seoul 
^ 9 University of Maribor, Maribor 
(Dated: May 30 2005) 

We report the results of a search for the flavor-changing neutral current process b — > d*y using 
a sample of 275 million B meson pairs accumulated by the Belle detector at KEKB. We find no 
significant signal for the exclusive decays B~ — > p~7, B° — > p°j, or B° — » uy. Assuming an 
isospin relation between the three modes, we set an upper limit for the combined branching fraction 
B(B — > (p, 01)7), of f.4 x 10 -6 at the 90% confidence level. This limit can be used to constrain the 
ratio of CKM matrix elements |Vtd/Vt s |. 

PACS numbers: 11.30.Hv, 13.40.Hq, 14.65.Fy, 14.40.Nd 



The b — > dq process, shown in Fig. Qfa), is a fla- 
vor changing neutral current transition that proceeds 
via loop diagrams in the Standard Model (SM). It is 
suppressed with respect to b — ► S7 by the Cabibbo- 
Kobayashi-Maskawa (CKM) factor \V td /Vts\ 2 ~ 0.04 
with a large uncertainty due to the lack of precise knowl- 
edge of Vtd- The exclusive modes B — > pj and B — > luj, 
which are presumably the easiest modes to search for, 
have not yet been observed 0. Calculations based on 
the measured rate for the b — ► S7 process B — > K*j 
that include IV^/l^sl 2 suppression, corrections due to 
form factors, St/ (3) breaking effects, and, for B~ decay, 
the additional annihilation diagram of Fig. ^b), pre- 
dict branching fractions in the range (0.9-2.7) x 10~ 6 
0, 0. Measurement of these exclusive branching frac- 
tions would improve the constraints on |Vtd/Vt s | in the 
context of the SM, and provide sensitivity to physics be- 
yond the SM that is complementary to that from b — > S7. 



(a) loop diagram (b) annihilation diagram 




FIG. I: (a) Loop diagram for b — » 1/7 and (b) annihilation 
diagram for B~ — » p~"/ only. 

In this paper, we report the results of a search for the 
b — > c?7 process using a sample of (275 ± 3) million B 
meson pairs accumulated at the T(4S) resonance. The 
data are produced in e + e~ annihilation at the KEKB 
energy-asymmetric (3.5 on 8 GeV) collider [4j and col- 



lected with the Belle detector 5]. The Belle detec- 
tor includes a silicon vertex detector (SVD), a central 
drift chamber (CDC), aerogel threshold Cherenkov coun- 
ters (ACC), timc-of-flight (TOF) scintillation counters, 
and an electromagnetic calorimeter (ECL) comprised of 
CsI(Tl) crystals located inside a 1.5 T superconduct- 
ing solenoid coil. An iron flux-return located outside of 
the coil is instrumented to identify muons (KLM). The 
dataset consists of two subsets: the first 152 million B 
meson pairs were collected with a 2.0 cm radius beampipe 
and a 3-layer SVD, and the remaining 123 million B me- 
son pairs with a 1.5 cm radius beampipe, a 4-layer SVD 
and a small-cell inner drift chamber |fj. 

We reconstruct the following final states: B~ — > p~7, 
B° — > p°7, and B° — » wy. Charge conjugate modes 
are implicitly included throughout this paper. We also 
reconstruct control samples of B~ — > K*~^ and B° — > 
K*°j decays. The following decay chains are used to 
reconstruct the intermediate states: p~ — > tt^tt , p° — > 
ir+ir-, uj -> tt+tt-tt , K*~ if^Tr , K*° -> R- tt+ , and 

7T° — * 77. 

Photon candidates are reconstructed from isolated 
clusters in the ECL that have no corresponding charged 
track, and a shower shape that is consistent with that of a 
photon. A photon in the barrel region of the ECL (33° < 
9 1 < 128° in the laboratory frame) with a center-of-mass 
(cm.) energy in the range 1.8 GeV < E* < 3.4 GeV is 
selected as the primary photon candidate. To suppress 
backgrounds from ir° — > 77 and r/ —> 77 decays, we veto 
the event if the reconstructed mass of the primary photon 
and any other photon of energy 30 (200) MeV or more is 
within a ±3<r window, i.e., ±18 (±32) MeV/c 2 , around 
the 7T° (77) mass. This set of criteria is referred to as the 
7T° /r) veto. For the primary photon, we sum the energy 
deposited in arrays of 3 x 3 and 5x5 ECL cells around 
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the maximum energy cell; if their ratio is less than 0.95, 
the event is vetoed. 

Neutral pions are formed from photon pairs with in- 
variant masses within ±10 (16) MeV/c 2 of the nominal 
7T° mass, corresponding to a ~2<r (^3er) window for the 
p~7 and K*~j (wy) modes. The photon momenta are 
then recalculated with a 7r° mass constraint. We require 
the energy of each photon to be greater than 30 MeV. Wc 
also require the cm. momentum of the 7r° to be greater 
than 0.5 GeV/c for the p~j and K*~j modes. 

Charged pions and kaons are reconstructed as tracks 
in the CDC and SVD. Each track is required to have a 
transverse momentum greater than 100 MeV/c and dis- 
tance of closest approach to the interaction point of less 
than 2 cm in radius and ±5 cm along the z-axis, which is 
aligned opposite to the positron beam. We do not use the 
track to form the signal candidate if, when combined with 
any other track, it forms a Kg candidate with an invari- 
ant mass within ±10 MeV/c 2 of the nominal Kg mass 
and a displaced vertex that is consistent with a K s . We 
determine the pion (L v ) and kaon (Lk) likelihoods from 
ACC, CDC and TOF information, and form a likelihood 
ratio L-xjiL-K + Lk) to separate pions from kaons. We 
select pions using criteria which have an efficiency of 89% 
and for which ~10% of kaons are misidentified as pions. 
For the <x>7 mode, we relax the requirement to select 94% 
of pions. (In the K*j modes, we select kaons with an ef- 
ficiency of 76-80%.) In addition, we remove tracks that 
are consistent with being electrons or muons. 

Invariant masses for the p and u> candidates are 
required to be within windows of ±150 MeV/c 2 and 
±30 MeV/c 2 , respectively, around their nominal masses. 

B candidates are reconstructed by combining a p 
or id candidate and the primary photon using two 
variables: the beam-energy constrained mass M^ c = 

\A-^bcam/ c2 ) 2 — \p*b/ c \ 2 an d tne energy difference AE = 
E* B — -Ebcami p* B and E B are the cm. momentum and 
energy of the B candidate, and E^ oam is the cm. 
beam energy. The magnitude of the photon momen- 
tum is replaced by {E£ eam — E*, u )/c when the mo- 
mentum p* B is calculated. To optimize the event selec- 
tion, we study Monte Carlo (MC) events in the region 
-0.10 GeV < AE < 0.08 GeV and 5.273 GeV/c 2 < 
Mbc < 5.285 GeV/c 2 : we choose selection criteria to 
maximize Ng / V Nb , where Ns is the MC signal yield 
assuming the SM branching fractions in Ref. |2| , and Nb 
is the expected background yield. 

The dominant background arises from continuum 
events (e + e~ — > 99(7)), where the accidental combina- 
tion of a p or u> candidate with a photon forms a B can- 
didate. We suppress this background using the following 
quantities: (1) J 7 , a Fisher discriminant Ml constructed 
from 16 modified Fox- Wolfram moments [a, E| and the 
scalar sum of the transverse momentum of all charged 
tracks and photons. (2) cos 9 B , where 9* B is the cm. polar 
angle of the B candidate direction: true B mesons follow 



a 1 — cos 2 0* B distribution, while candidates in the contin- 
uum background are uniformly distributed. (3) Az, the 
separation between the decay vertex of the candidate B 
meson and the origin of the remaining tracks in the event 
along the z-axis. The two vertices are reconstructed in 
about 85% of the p°j and 077 events. Discrimination 
is provided due to the displacement of the signal B de- 
cay vertex from the other B, as tracks from continuum 
events typically have a common vertex. (4) A tagging 
quality r which indicates the level of confidence in the 
-B-flavor determination described in Ref. 10]. The algo- 
rithm provides additional discrimination between signal 
and continuum background where no true B meson is 
present. 

For each of the quantities J 7 , cos0 B and Az, we con- 
struct one-dimensional likelihood distributions for sig- 
nal and continuum events. Signal distributions are 
modeled with an asymmetric Gaussian function for J 7 , 
|(ao — a 2 cos 2 9* B ) for cos 9* B , and an exponential con- 
volved with a Gaussian resolution function for Az. Con- 
tinuum background distributions are modeled with an 
asymmetric Gaussian for J 7 , (b — b 2 cos 2 9* B ) for cos8* B , 
and a sum of three Gaussians with a common mean for 
Az; the coefficients eto, a 2l bo are close to unity while 
the coefficient b 2 is close to zero. The Az likelihood 
distribution for continuum is determined from data in 
the sideband region 5.20 GeV/c 2 < M hc < 5.24 GeV/c 2 , 
\AE\ < 0.3 GeV; we treat the two subsets of the data 
separately, as the vertex resolution is improved in the 
subset with the 4-layer SVD. The T and cos#g likeli- 
hood distributions are determined from MC samples. 

We form product likelihoods C s and C c for signal and 
continuum background, respectively, from the likelihood 
distributions for J 7 , cos9 B and (where available) Az, 
since these quantities are independent. In the plane 
(r, 7V) defined by the tagging quality r and the likeli- 
hood ratio 1Z = C S /(C S + C c ), signal tends to populate 
the edges at r — 1 and TZ — 1, while continuum tends to 
populate the edges at r — and 1Z = 0. We select the 
events in a signal enriched region defined by 1Z > 1Z± for 
r > rt, and 1Z > 1 — a(l + r) for r 2 < r < r±, where 
the parameters r±, r 2l TZi and a are mode- and subset- 
dependent and are determined so that Ng/^Ng + Nb 
is maximized (we use this quantity instead of Ns/^/Nb 
because of the limited size of the MC simulation sam- 
ple used in this procedure). The values are r\ = 0.85, 
r 2 = 0.01, Tlx = 0.8 and a = 0.025 for the first subset of 
the /9°7 mode; similar values are used for the other sub- 
sets. We define the rest of the area as the background 
enriched region. 

We consider the following B decay backgrounds: B — > 
K*j, other B — > processes, B — > pir° and ujtt , 

B — * prj and ujij, B~ — > p~ p , other charmless B decays, 
and the b — > c transition. We find the b — > c back- 
ground to be negligible. The B — > K*j background may 
mimic the signal decay B — > p"f if the kaon from K* is 
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misidentified as a pion. To suppress B — > if* 7 events 
we calculate Mk-k, where the kaon mass is assigned to 
one of the charged pion candidates, and reject the can- 
didate if M K7T < 0.96 (0.92) GeV/c 2 forthe p°7 (p-j) 
mode. The decay chain B a -> if*°7, if* -► if^Tr , 
if <y — ► 7r + 7r _ has a small contribution to B° — > cl>7 due 
to the tail of the if* Breit-Wigner lineshape. In addi- 
tion, B — * K*j and other B — > A s 7 decays contribute to 
the background when the p and w candidates are selected 
from a random combination of particles. 

Charmless decays B — > pir°, lutt , prj, and uirj may 
mimic the signal if a photon from 7r° or 77 — > 77 decay is 
soft and undetected by the 7r°/77 veto. To suppress this 
background, we reject the candidate if |cos#hei| > 0.8 
(0.6) for the p°7 and W7 (p _ 7) modes, where the helicity 
angle #hel is the angle between the 7r + and B momentum 
vectors in the p rest frame, or the angle between the 
normal to the u> decay plane and the B momentum vector 
in the u> rest frame. The decay B~ — > p~ p , p~ — > 7r~7r° 
also contributes to the B a — > p°7 mode when a pion from 
the p~ decay and a photon from the tt° decay are both 
soft and undetected. Other charmless decays have very 
small contributions and are considered as an additional 
background component when we extract the signal yield. 

The reconstruction efficiency for each mode is defined 
as the fraction of the signal yield remaining after all se- 
lection criteria, where the signal yield is determined from 
an extended unbinned maximum likelihood fit to the MC 
sample. Each signal distribution is modeled as the prod- 
uct of a Gaussian in M^ c and a Crystal Ball lineshape [ll| 
in AE to reproduce the asymmetric ECL energy response 
for AE. The background component is modeled as the 
product of a linear function for AE and an ARGUS func- 
tion ^3 f° r Mhc- The efficiencies are listed in Table [I] 
The systematic error on the efficiency is the quadratic 
sum of the following contributions, estimated using con- 
trol samples: the photon detection efficiency (2.2%) as 
measured in radiative Bhabha events; charged tracking 
efficiency (1.0% per track), from partially reconstructed 
D*+ -> D°ir + , D° -> K a s TT+ir-, K° s -> 7r+(7r-); charged 
pion identification (1.0% per pion) from D* + — ► D°ir + , 
D° — > if ~7r + ; neutral pion detection (4.6-7.3%) from ij 
decays to 77, 7r + 7r~7r° and Sir ; IZ-r and 7r°/ry veto re- 
quirements (5.4%) from B~ — > D°ir~ , D° — > K~ir + ; and 
uncertainty due to limited MC statistics (0.9-1.5%). 

We perform an unbinned maximum likelihood fit to 
the data in the (M bc , AE) region bounded by \AE\ < 
0.3 GeV and Mb c > 5.2 GeV/c 2 , simultaneously for the 
three signal modes and the two B — ► K*-f modes; the lat- 
ter are used to determine the residual K*j backgrounds 
in the signal modes. We fit the two data subsets simul- 
taneously, so that in total ten distributions are included 
in the fit. We define the combined branching fraction 
B(B — ► (p, 10)7) = B(B~ — > p~7), assuming the isospin 
relation [lj B(B~ -> p'j) = 2^B(B° -> p° 7 ) = 



2Ls±B(B° -> w 7 ); we use ^± = 1.086 ± 0.017 0. 

We also assume B{B -> K*j) = B(B~ -> K*~j) = 
^B(B -> if *° 7 ). 

We describe the events in the fit region using a sum 
of functions for the signal, continuum, if *7, and other 
background hypotheses. Signal parameters for M^ c and 
AE are calibrated using B~ — > D°ir~ and £? — > If* 7 
samples, respectively. We use the combined branching 
fractions -B — * (p, w)7 and B — > if* 7 as free parameters. 
The continuum background is modeled as the product of 
a linear function in AE whose slope is allowed to float, 
and an ARGUS function in Mb c whose parameters are 
fixed from a comparison between data and MC in the 
background enriched region and MC in the signal en- 
riched region. The continuum contribution in the data is 
allowed to float. The size of the if* 7 background compo- 
nent in each (p, lo)j channel is constrained using the fit to 
the if *7 events. The levels of the other backgrounds are 
fixed using known branching fractions or upper limits. 

Results of the fit are shown in Fig. El and listed in Ta- 
ble 0] The significance of the fit, •*/— 2 ln(£ / £max) , is 
found to be 1.2tr; where C max (£0) is the maximum like- 
lihood from the fit when the signal branching fraction 
is floated (set to zero). To include the effect of possible 
systematic error in this calculation, we change each pa- 
rameter by ±lcr in the direction that gives the smaller 
resulting significance. The systematic error in the signal 
yield is estimated by varying each of the fixed param- 
eters by lcr, and then taking the quadratic sum of the 
deviations in the signal yield from the nominal value. 
The combined branching fraction is B(B — > (p,uj)^) = 
(0.72 ± 0.28) x 10~ 6 , where the first and second er- 
rors are statistical and systematic, respectively. 

Since the significance is small, we quote a 90% 
confidence level upper limit Bgo using the formula 
/ 90 C{x)dx = 0.9 f °° C(x)dx, where C(x) is the likeli- 
hood function with the branching fraction fixed at x. The 
systematic error is taken into account assuming a Gaus- 
sian distribution. We find B(B — > (p,uj)-f) < 1.4 x 10~ 6 
at the 90% confidence level. We also perform individual 
fits to the three signal modes; the corresponding upper 
limits are listed in Tabic [I] 

A similar fit is performed using the ratio of branch- 
ing fractions B(B — > (p,uj)^)/B(B — > if* 7) instead of 
B(B — > (p, u))y), so that the systematic error partially 
cancels. We find B(B -> ( P ,lu)j)/B{B -► if* 7 ) < 0.035 
at the 90% confidence level. One can use this result to 
constrain Vtd'- for example, using the prescription given 




AR] where the form factor ratio ( — 0.85 ± 0.10 
and S'C/(3)-breaking effect AR = 0.1 ± 0.1, we obtain 
\V t d/V ts \ < 0.22 at the 90% confidence level. This limit 
is slightly smaller than that resulting from searches for 
B s mixing 14]. 



In conclusion, we search for the b — > process using 
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a simultaneous fit to the B — > (p, u)j and B — > if* 7 
modes. The upper limit we obtain is within the range 
of SM predictions 0, and can be used to constrain 
\V td /V ts \. 
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FIG. 2: Projections of the simultaneous fit results to Mt, c 
(in the region -0.10 GeV < AE < 0.08 GeV) and AE (in 
the region 5.273 GeV/c 2 < M bc < 5.285 GeV/c 2 ) for the 
individual modes and their sum. Lines represent the total 
fit result (solid), signal (dashed), continuum (dotted), and B 
decay background (dot-dashed) components. 
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TABLE I: Yield, significance (including systematic error), efficiency and fitted branching fraction (central value and 90% 
confidence level upper limit) for each mode and for B — ► (p,u))~f. 

Branching Fraction 



Mode 


Yield 


Significance 


Efficiency 


Central Value 


Upper Limit 


B~ -> p~-y 
W p° 7 
B° -> twy 
Combined 


18.7 +1 °i 
1.91" 

o-9±|:l 


1.7a 
0.3cr 

0.1(7 
1.2(7 


(5.5 ± 0.4)% 
(3.9 ± 0.3)% 
(3.9 ± 0.4)% 


(1.24 t° ;ll ±0.26) x 10~ 6 
(0-18 toils ±0.10) x 10~ 6 
(0.08 ±g;|? ±0.19) x 10~ 6 
(0.72 ±° *l ± 0.28) x 10~ 6 


2.2 x 10 -6 
0.8 x 10~ 6 
0.8 x 10~ 6 
1.4 x 10~ 6 



